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Abstract 

The stopping power of protons and deuterons in low energy collisions with helium gas targets 
is investigated with the numerical solution of the time-dependent Schrodinger coupled-channels 
equations using molecular orbital wavefunctions. It is shown that at low projectile energies the 
energy loss is mainly due to nuclear stopping, charge exchange, and the excitation of low energy 
levels in the target. The second and third mechanisms, called electronic stopping, dominate for 
E/afe > 200 eV. At lower energies it is also shown that a threshold effect is responsible for a quick 
drop of the energy loss. This investigation sheds more light on the long standing electron screening 
problem in fusion reactions of astrophysical interest. 
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Nuclear fusion reactions proceed in stars at extremely low energies, e.g. of the order of 
10 keV in our sun . At such low energies it is extremely difficult to measure the cross 
sections for charged particles at laboratory conditions due to the large Coulomb barrier. 
Moreover, laboratory measurements of low energy fusion reactions are strongly influenced 
by the presence of the atomic electrons. One has observed experimentally a large discrepancy 
between the experimental data and the best models available to treat the screening effect 
by the electrons in the target nuclei j^]. The screening effect arises because as the projectile 
nucleus penetrates the electronic cloud of the target the electrons become more bound and 
the projectile energy increases by energy conservation. Since fusion cross sections increase 
strongly with the projectile's energy, this tiny amount of energy gain (of order of 10-100 
eV) leads to a large effect on the measured cross sections. However, in order to explain 
the experimental data it is necessary an extra-amount of energy about twice the expected 
theoretical value [3|. 

In order to extract the fusion cross sections from experiment one needs to correct for 
the energy loss in the target to assign the correct projectile energy value for the reaction. 
The authors in refs. Q and Q have shown that a possible solution to the long standing 
discrepancy between theory and experiment for the reaction 3 He(d, p) 4 He could be obtained 
if the projectile energy loss by electronic excitations and charge exchange with the target 
atoms would be smaller than previously assumed in the experimental data analysis. There 
have been indeed a few experiments in which evidences of smaller than expected electronic 
stopping power were reported (see, e.g. ref. Other reactions of astrophysical interest 

(e.g., those listed in by Rolfs and collaborators 0, 3|) should also be corrected for this 
effect. Whereas at higher energies the stopping is mainly due to the ionization of the target 
electrons, at the astrophysical energies it is mainly due to excitations of the lowest levels, 
charge-exchange between the target and the projectile, and the nuclear stopping power. 

In this work I address the problem of the stopping of very low energy ions in matter. I 
consider the systems p+ 4 He and d+ 3 He, for which there are experimental data available. A 
previous work studied the energy loss of protons on hydrogen gas targets and showed that 
the stopping at very low proton energies is indeed smaller than what would be expected from 
extrapolations based on the Andersen and Ziegler tables 10]. The case of helium targets is 
more complicated due to the electron-electron interaction. 

The present approach is based on the solution of the time-dependent Schrodinger equation 
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for the electron in a dynamical two-center field. The transition from the separated atoms 
(H + + He) and the united atom (Li + ) is obtained in the adiabatic approximation, i.e. by 
assuming that the electronic motion is fast compared to the nuclear separation motion so that 
the molecular orbitals (MO) are those for the distance R(t) between the nuclei. The atomic 
wavefunctions, A, = J2j c j^ ( t ) f at i are constructed as a linear combination of Slater-type 
orbitals (STO) [ll| of the form (j)^ lat = iVY n_1 exp (— £r) Y\ m (r), where the Slater coefficients 
n and ( are chosen to best approximate the exact atomic wavefunctions (see, e.g. ref. 
The molecular orbital wavefunctions for the p + He system, are obtained with the 0^'s chosen 
so that half of the STO's are centered on the proton (A) and the other half are centered 
on the helium nucleus (B). The total wavefunction for the two-electron system is finally 
written as a Slater determinant of the molecular orbital wavefunctions, 
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where a, [3 denote the spin state of the electron. Configuration-interaction with double 
excitation configurations were included in the calculation, with the coefficients n and the 
Slater parameters ( chosen in a variational method to obtain the lowest energy states of the 
system. 

Using these conditions and variation method, one obtains the following Hatree-Fock equa- 
tion: F<g>C = 0<g>C<g>E, where F is the "Fock" matrix 
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in which "occ" refers to the occupied molecular orbital 
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are the two-electron integrals. The C matrix is the coefficient matrix and O is the overlap 
matrix (0 M (1) \<p v (1)) . E is a diagonal matrix with each diagonal element corresponding to 
the energy of the associated molecular orbital. Solving the Hartree-Fock equations one 
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obtains the coefficients Cy which give the proper linear combination of atomic orbitals to 
form the molecular orbital. The energy of the molecular orbitals are then given by 



E(R) = P^ H ^ + \Ys 





1 








r\2 



Act 



1 








r\2 





(5) 



Table I shows the states involved in the calculation where it is shown how the states 
in the separated hydrogen and helium atoms become molecular states in the united atom 
system. For large distances between the nuclei, R > 15 a.u. (1 a.u. of length = 0.53 A) the 
energy levels for the Is, 2s, and 2p states of H and He are reproduced to within 2% and 4% 
of the spectroscopic data, respectively. The energies of these states are shown in figure [T] as 
a function of the internuclear distance R. 



Separated atom 


United atom 


H+ + He(ls 2 ) 


OS 


H(ls) +He + (ls) 


IS 


H+(ls) + He(ls2s) 


2S 


H(n = 2) + He + (ls) 


in 


H(n = 2) + He + (ls) 


3S 


H(n = 2) + He + (ls) 


4S 


H+ + He(lslp) 


5S 


H+ + He(lslp) 


2n 



Table 1 - Lowest states in the p+He mo 
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At very low proton energies (Ep < 10 keV) it is fair to assume that only the low-lying 
states are involved in the electronic dynamics. Only for bombarding energies larger than 
25 keV the proton velocity will be comparable to the electron velocity, v e ~ ac. Thus, 
the evolution of the system is almost adiabatic at E p < 10 keV. Also shown in figure ^ 
(inset) are the intersection points of the states with same symmetry. In a fast collision these 
states would cross (diabatic collisions), whereas in collisions at very low energies (adiabatic 
collisions) they obey the von Neumann- Wigner non-crossing rule. 

In the dynamical case the full time-dependent wavefunction for the system can be ex- 
panded in terms of two-center states, ip n (ri, r 2 , t), given by eq. HJ with expansion coefficients 
a n (t). It is further assumed that the nuclei follow a classical straight-line trajectory deter- 
mined by an impact parameter b, so that the time dependence of the molecular wavefunctions 
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FIG. 1: Adiabatic energies (1 a.u. of energy = 27.2 eV, 1 a.u. of length = 0.53 A) for the electronic 
orbitals for the (H-He) + system as a function of the internuclear separation. As the atoms approach 
each other slowly curves of same symmetry repel each other. A transition between states s and s' 
can occur in a slow collision. In a fast collision a diabatic transition, with the states crossing each 
other, will occur. This is shown in the inset. 




is determined by the condition R = \JW- + vH 2 , where v is the collision velocity. The dy- 
namical evolution of the H+He system is calculated using the same approach as described 
in ref . 9| . We solve the set of linear coupled equations 

dA 

tS-^=M-A (6) 

dt v ; 

where the column matrix A represents the time-dependent expansion coefficients, S is the 
overlap matrix with elements S^- = and M is the coupling matrix with elements 

Mij = (ipi\H e i — id / dt\ , where H e i is the electronic Hamiltonian. The solutions are 
obtained starting from initial internuclear distance of 15 a.u. for the incoming trajectory 
and stopped at the same value for the outgoing trajectory. The probability for the capture 
in the proton is obtained by a projection of the final wavefunction into the wavefunctions of 
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the Is, 2s and and 2p states of the hydrogen atom. 
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Resonant charge-exchange in atomic collisions was first observed by Everhart and col- 



laborators 



12j. In these experiments it was determined that the exchange probability in 



homonuclear atomic collisions oscillates with the incoming energy for a collision with a 
given impact parameter, or scattering angle. This was interpreted [13| as due to transitions 
between degenerate states of the system at large internuclear separation distance. In the 
simplest situation of a p + H the degenerate states are the symmetric and antisymmetric 
states obtained from the linear combination of the (H + H) and (HH + ) wavef unctions. This 



effect was studied in ref. jl4| . where a relation between the damping of the oscillatory 
behavior of the exchange probabilities and the Landau-Zener effect was established. The 
p + H collisions at small energies was recently studied in ref. [IJ and the oscillatory effect 
was shown to be related to the Sommerfeld quantization rule for the integral from t = — oo 
to t = — oo of the energy difference between the symmetric and antisymmetric state. The 
electron tunnels back and forth between the projectile and the target during the ingoing and 
the outgoing parts of the trajectory. When the interaction time is an exact multiple of the 
oscillation time, a minimum in the exchange probability occurs. 

A similar situation occurs for p + He collisions, as shown in figure 121 for the electron 
capture probability by the proton at 10 keV bombarding energy. These oscillations are 
due to the electron exchange between the ground state of the hydrogen and the first excited 
state in He (ls2s). But, in contrast to the H + H system, the oscillations are strongly damped. 
Following the work of Lichten Q| we interpret this damping effect as due to the interference 
between the participant states and a band of states of average energy (E a ) and width 2T, 
as seen in figure ^ The important regions where the diabatic level cross occurs is shown 
in figure Q inside the encircled areas. The damping mechanism is best understood using 
the Landau-Zener theory for level crossing. At the crossing there is a particular probability 
(1 — P) of an adiabatic transition where P is given by the Landau-Zener formula 

Pexch = 6XP y(d/dR) (E 8 -E s ,)_ (8) 
where v is the collision velocity and H ss / is the off-diagonal matrix element connecting states 
s and sf. The oscillatory behavior shown in figure El is due to the many level transitions 
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FIG. 2: Probability of charge exchange in the collision p+ 4 He showing the resonant behavior as a 
function of the impact parameter and for proton energy E p = 10 KeV. 

at the crossing, each time governed by the probability of eq. |H1 The interference with the 
neighboring states introduces a damping in the charge exchange probability, i.e. 



where (E a ) ~ 1. a.u. is the average separation energy between the OS level and the bunch 
of higher energy levels shown in figure ^ The exponential damping factor agrees with the 
numerical calculations if one uses T ~ 5 eV, which agrees with the energy interval of the 
band of states shown in figure H 

The total cross section for charge exchange is calculated from 



The numerical results for the p + He system is shown in figure 01 as a function of the proton 
energy. The solid line is the result of using the coupled-channels equations H3 and [7| for 
the exchange probability. The experimental data are from ref. [15]. We observe that the 
calculation reproduces the trend of the experimental data. But, whereas the maximum of 
the cross section at E p ~ 20 keV is rather well described, the calculations underestimate the 
cross sections at smaller energies. The low energy slope of the cross section is nonetheless 
well reproduced. At energies higher than the Bragg peak (E p > 20 keV) the numerical results 
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FIG. 3: Charge-exchange cross sections for the p + 4 He as a function of the proton energy. The 
solid line was obtained by solving the coupled-channels equations El and using for the exchange 
probability. The experimental data are from ref. 



should not be trusted as the adiabatic approximation for the molecular orbitals and also the 
inclusion of only the lowest energy levels are not adequate (continuum states should also be 
included). For E p — > 0, the charge exchange cross section must go to zero since the higher 
binding of the electrons in He prevents the capture by the incident proton in an extreme 
adiabatic collision. This feature is correctly reproduced by the numerical calculations. 

In figure HI we show the stopping cross section of the proton. The experimental data are 
from refs. Q]- The stopping cross section is denned as S = J2i , where AEi is 

the energy loss of the projectile in a process denoted by i. The stopping power, Sp = dE/dx, 
the energy loss per unit length of the target material, is related to the stopping cross section 
by S = Sp/N, where N is the atomic density of the material. In the charge exchange 
mechanism one of the electrons in He is transferred to incoming proton and the energy loss 
by the proton is given by AE = m e v 2 /2, where v is the proton velocity. Assuming that 
there are few free electrons in the material (e.g. in the helium gas) only one more stopping 
mechanism at very low energies should be considered: the nuclear stopping power, S n . This 
is simply the elastic scattering of the projectile off the target nuclei. The projectile energy 
is partially transferred to the recoil energy of the target atom. The stopping cross section 
for this mechanism has been extensively discussed in ref. [l2|. The total stopping power is 
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FIG. 4: Stopping cross section for proton incident on gas 4 He targets, as a function of the proton 
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energy. The experimental data are from refs. 



given by S = S exc h + S n . In units of 10 5 eV cm 2 the nuclear stopping for the p+ 4 He system 
at E p < 30 keV is given by 

In (l + 1.1383e) 
(e + 0.0132Ls°- 21226 + 0.19 5 93£ - 5 ) ' U 

where So — 0.779 and e = 5.99E P , with the proton incident energy E p given in keV. 

The dashed line in figure H] gives the energy transfer by means of nuclear stopping, while 
the solid line shows the results for the electronic stopping mechanism, i.e. due to charge- 
exchange and excitation in the helium target. We see that the nuclear stopping dominates 
at the lowest energies, while the electronic stopping is larger for proton energies greater than 
200 eV. We do not consider the change of the charge state of the protons as they penetrate 
the target material. The exchange mechanism transforms the protons into H atoms. These 
again interact with the target atoms. They can loose their electron again by transferring it 
back to a bound state in the target. 

At very low energies the only possibility that the electron is captured by the proton is 
if there is a transition ls 2 ( 1 So) — ► ls2s( 3 S) in the helium target. Only in this case the 
energy of one of the electrons in helium roughly matches the electronic energy of the ground 
state in H. This resonant transfer effect is responsible for the large capture cross sections. 
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When this transition is not possible the electrons prefer to stay in the helium target, as 
the energy of the whole system is lowest in this case. Another possible mechanism for the 
stopping is the excitation of the helium atom by the transition ls 2 ( 1 S ) — > ls2s( 3 S). Thus, 
there must be a direct relationship between the energy transfer to the transition ls 2 ( 1 Sn) 
— > ls2s( 3 S) and the minimum projectile energy which enables electronic changes. Ref. 18[ 
reported for the first time this effect, named by threshold energy, which can be understood 
as follows. The momentum transfer in the projectile-target collision, Aq, is related to the 
energy transfer to the electrons by Aq = AE/v, where v is the projectile velocity. In order 
that this momentum transfer absorbed by the electron, induces an atomic transition, it is 
necessary that h 2 Aq 2 /2m e ~ AE. Solving these equations for the projectile energy one finds 

E thres _ UhLAE . (10) 

p 4m e v ' 

This is the threshold energy for atomic excitations and/or charge exchange. If the projectile 

energy is smaller than this value, no stopping should occur. The energy for transition 

ls 2 ( 1 So) — > ls2s( 3 S) in He is AE = 18.7 eV. Thus, for p + He collisions, the threshold 

energy is E t ^ ires ~ 9 keV. This roughly agrees with the numerical calculations presented in 

figure U] (solid curve). 

Figure El shows the e nerg y loss of deuterons in 3 He g function of deuteron energy. 

The data are from ref. |l8[. The solid curve is the numerical calculation for the electronic 
stopping power, while the dashed curve shows the nuclear stopping. For this system the 
coefficients in eq. 01 are Sq = 1.557 and e = 4.491.2^, respectively. As discussed in ref Q| 
the threshold deuteron energy in this reaction is of the order of 18 keV, which agrees with 
the estimate based on eq. El However, the numerical calculations based on the electronic 
stopping (solid curve of fig. EJ) indicate a lower threshold energy for this system. Nonetheless, 
the agreement with the experimental data is very good for E^ > 20 keV. The threshold effect 
is one more indication that the extrapolation S ~ v, based on the Andersen- Ziegler tables 
is not applicable to very low energies. 

The steep rise of the fusion cross sections at astrophysical energies amplifies all effects 
leading to a slight modification of the projectile energy [19]. The results presented in this 
article show that the stopping mechanism does not follow a universal pattern for all systems. 
The threshold effect reported in ref. [l^ is indeed responsible for a rapid decrease of the 
electronic stopping at low energies. It will occur whenever the charge-exchange mechanism 
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^IG. 5: Energy loss of deuterons in 3 He gas as a function of deuteron energy. Data are from ref. 
l3 | . The solid curve is the calculation for the electronic stopping power, while the dashed curve 
shows the nuclear stopping. 



and the excitation of the first electronic state in the target involve approximately the same 
energy. However, the drop of the electronic stopping is not as sharp as expected from the 
simple classical arguments given by eq. fTUl 

The experiments on astrophysical fusion reactions have shown that the screening effect 
is much larger than expected by theory. The solution to this problem might be indeed the 
smaller stopping power, due to a steeper slope at low energies induced, e.g. by the threshold 
mechanism. This calls for improved theoretical studies of the energy loss of ions at extremely 
low energies of and for their independent experimental verification. The present situation is 
highly disturbing because if we cannot explain the laboratory screening effect, most likely 
we cannot explain it in stellar environments. 
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